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Central issues in the electronic structure of underdoped cuprate superconductors are 
to clarify the shape of the Fermi surfaces and the origin of the pseudogap. On the 
basis of the model proposed by Kamimura and Suwa, which bears important features 
originating from the interplay of Jahn- Teller physics and Mott physics, the feature 
of Fermi surfaces in underdoped cuprates is the presence of Fermi pockets constructed 
from doped holes under the coexistence of a metallic state and a local antiferromagnetic 
order. Below T^, the holes on Fermi pockets form Cooper pairs with d-wave symmetry 
in the nodal region. In the antinodal region, there are no Fermi surfaces. In this study 
we calculate the energy distribution curves (EDCs) of angle-resolved photoemission 
spectroscopy (ARPES) below T^. It is shown that the feature of ARPES profiles of 
underdoped cuprates consists of a coherent peak in the nodal region and real transitions 
of photoexcited electrons from occupied states below the Fermi level to a free-electron 
state above the vacuum level in the antinodal region, where the latter transitions form a 
broad hump. From this feature, the origin of the two distinct gaps observed by ARPES 
is elucidated without introducing the concept of the pseudogap. Finally, a remark is 
made on the phase diagram of underdoped cuprates. 
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Undoped copper oxide (La2Cu04) is an antiferromagnetic Mott insulator, in which 
an electron correlation plays an important role.^^ Thus, we may say that undoped 
cuprates are governed by Mott physics. In 1986, Bednorz and Miiller discovered high- 
temperature superconductivity in copper oxides by doping hole carriers into La2Cu04.^^ 
Their motivation was the consideration that higher Tc could be achieved for copper oxide 
materials by combining Jahn- Teller (JT) active Cu ions with the structural complex- 
ity of layer-type perovskite oxides. To investigate the mechanism of high-temperature 
superconductivity, it is assumed in most models that doped holes itinerate through 
orbitals extending over a Cu02 plane in systems consisting of CuOe octahedrons elon- 
gated by the JT effect. These models are called the "single-component theory", because 
the orbitals of hole carriers extend only over a Cu02 plane. 

In 1989, Kamimura and coworkers showed by first-principles calculations that the 
apical oxygen in CuOg octahedrons tends to approach Cu^+ ions when Sr^+ ions are 
substituted for La^"*" ions in La2Cu04 in order to gain the attractive electrostatic energy 
in ionic crystals such as cuprates.^' As a result, CuOg elongated by the JT effect 
shrinks with hole doping. This deformation against the JT distortion is called the"anti- 
Jahn- Teller effect".^) By this effect, the energy separation between the two kinds of 
orbital states, which have been spht originally by the JT effect, becomes smaller with 
hole carrier doping. These two states are the aig antibonding orbital state |a*g) and big 
bonding orbital state |big). The a*g antibonding orbital state is constructed by a Cu dz"^ 
orbital and the six surrounding oxygen p orbitals including apical O p^-orbitals; the |big) 
orbital is constructed by four in-plane O Po- orbitals with a small Cu d,-t.2_y2 component 
parallel to the Cu02 plane. The spatial extensions of the ja^g) and |big) orbitals, which 
are perpendicular and parallel to the Cu02 plane, respectively, are shown in Fig. 1. 

By taking account of the anti- Jahn- Teller effect, Kamimura and Suwa reported that 
one must consider these two kinds of orbital states equally in forming the metallic state 
of cuprates; they constructed a metallic state coexisting with the local antiferromagnetic 
(AF) order. ^) This model is called the "Kamimura-Suwa (K-S) model". Since the anti- 
Jahn- Teller effect is a central issue of Jahn- Teller physics, we may say that the K-S 
model bears important features originating from the interplay of Jahn- Teller physics 
and Mott physics. Since these two kinds of orbitals extend not only over the Cu02 plane 
but also along the direction perpendicular to it, the K-S model represents a prototype 
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Fig. 1. (Color online) Spatial extension of a^g antibonding orbital ja^g) and big bonding orbital 
|big}. a, a|g antibonding orbital |a|g). b, big bonding orbital |big). 

of a "two-component theory" , in contrast to the single- component theory. 

On the basis of the K-S model, Kamimura and Ushio have calculated Fermi surfaces 
in underdoped La2-.a;Sr^Cu04 (LSCO),^'^^ and have shown that the coexistence of a 
metallic state and a local AF order results in the Fermi pockets constructed from doped 
holes in the nodal region. The appearance of Fermi pockets and small Fermi surfaces 
in cuprates has recently been reported by various experimental groups. 

In this study, on the basis of the K-S model, we calculate the energy distribution 
curves (EDCs) of the angle-resolved photoemission spectroscopy (ARPES) profiles of 
cuprates below T^, and we show that the feature of the calculated ARPES profiles 
consists of a coherent peak due to the superconducting density of states in the nodal 
region and the real transitions of electrons from the occupied states below the Fermi level 
to a free-electron state above the vacuum level in the antinodal region. In particular, 
we show that the latter transitions form a broad hump in ARPES EDCs in underdoped 
cuprates. 

Concerning the ARPES experiments in underdoped cuprates, Tanaka and coworkers 
reported very interesting gap features in their observation of ARPES spectra. Their 
result exhibits a coherent peak in the nodal region and a broad hump in the antinodal 
region in underdoped Bi2212 samples below T^,.^^^ From the quantitative agreement 
between the theory and the experiment, we conclude that the observed broad hump 
corresponds to the photoelectron excitations from the occupied states below the Fermi 
level to the free-electron state above the vacuum level. In this context, it is concluded 
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Finally, in connection with the finite size of the spin-correlation length in a metallic 

state, we discuss the finite size effect of a metallic state on the spin-electronic structures 

of underdoped cuprates, and a new explanation for the phase diagram for underdoped 

cuprates is proposed. 

The organization of the present paper is as follows: At the beginning of §2 we first 

summarize the essential features of the K-S model, which bears important features 

originating from the interplay of JT physics and Mott physics. In §3, on the basis of the 

many-body effects including energy bands obtained from the K-S model, we predict the 

key features of ARPES EDCs and clarify the origin of the two-gap scenario proposed 

from the experimental results of Tanaka et al}^^ In §4, we discuss the finite size effects 

on the Fermi surfaces in cuprates. In connection with the finite size effects, we discuss 

the possibility of the spatially inhomogeneous distribution of Fermi pocket states and 

large Fermi surface states. Taking account of the finite size effect, we propose a new 

interpretation for the phase diagram of underdoped cuprates in §5. We devote §6 to the 

conclusion and concluding remarks. 

2. On the K-S Model 

In this section, we summarize the main features of the K-S model, emphasizing its 
important roles in underdoped cuprates due to the interplay of JT physics and Mott 
physics. 

2. 1 Key features of the K-S model 

The key features of the K-S model are explained in a heuristic way using the picture 
of a two-story house model shown in Fig. 2. In this figure, the first story of a Cu 
house is occupied by Cu localized spins, which form the AF order in the spin-correlated 
region by the superexchange interaction. The second story in the Cu house consists of 
two fioors due to the anti-JT effect; the lower a^g floor and the upper big floor. The 
second stories of neighboring Cu houses are connected by oxygen rooms, reflecting the 
hybridization of Cu d and O p orbitals. In the second story, a hole carrier with an 
up-spin enters the a*g floor of the left-hand Cu house owing to Hund's coupling with a 
Cu localized up-spin in the first story (Hund's coupling triplet),^''' as shown in the 
leftmost column of the figure. By the transfer interaction marked by a long arrow in the 
figure, the hole is transferred into the big fioor in the neighboring Cu house (second from 
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Fig. 2. (Color online) Explanation of K-S model using picture of two-story house model. 

the left) through the oxygen rooms, where a hole with up-spin forms a spin-singlet state 
with a localized down-spin in the second Cu house from the left (Zhang- Rice singlet). ^^-^ 
The key feature of the K-S model is that the hole carriers in the underdoped regime 
of LSCO form a metallic state by taking the Hund coupling triplet and the Zhang- 
Rice singlet alternately in the presence of a local AF order without destroying the AF 
order, as shown in the figure. From Fig. 2, one may understand that the characteristic 
feature of the K-S model is the coexistence of the AF order and a normal, metallic 
(or a superconducting) state in the underdoped regime. This feature of the K-S model 
(two-component theory) is different from that of the single-component theory. 

As seen in Fig. 2, the wave functions of a hole carrier with up and down-spins have 
the following phase relation: 

%(^) =exp(^^-a)%(r-). (1) 

Kamimura et al. have shown that this unique phase relation leads to the d-wave super- 
conductivity. 

2.2 Effective Hamiltonian for the K-S model 

The following effective Hamiltonian is introduced to describe the K-S model follow- 
ing Kamimura and Suwa^-* (see also ref. 7). It consists of four parts: the one-electron 
Hamiltonian i^smg for the a*g and big orbital states, the transfer interaction between 
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dx2_y2 localized spins Hp^p, and the exchange interactions between the spins of dopant 

holes and dx2_y2 localized holes within the same CuOe octahedron iJex- Thus, we have 

H — -fusing + + i^AF + -f^ex 

i,m,a 

+ tmn {pima^jna + h.C.^ 

{i,j>,m,n,(T 

+ J ^ ^ Si ■ Sj + ^ ^ Si^rn ■ 'S'i , (2) 

{i,j) i,m 

where (m = a^g or big) represents the one-electron energy of the aj^g and big orbital 
states, Cj^^ and Ci^a are the creation and annihilation operators of a dopant hole 
with spin a in the ith CuOg octahedron, respectively, t^n is the transfer integral of a 
dopant hole between the m-type and n-type orbitals of neighboring CuOe octahedrons, 

— * — * 

J is the superexchange interaction between spins Si and Sj of d^^^yi localized holes 
in the b*g orbital in the nearest- neighbor Cu sites i and j (J > for AF interaction), 
and K„i is the exchange integral for the exchange interactions between the spin of a 
dopant hole, Sim, and the d^-i-y-i localized spin Si in the ith CuOg octahedron. There 
are two exchange constants, i.e. i^Ta*^ and -ft^big, for the Hund coupling triplet and the 
Zhang- Rice singlet, respectively, where i^a*^ < and -f^big > 0. The appearance of the 
two kinds of exchange interactions in the fourth term is due to the interplay of Mott 
physics and JT physics. This is the key feature of the K-S model. 

The electron-electron interactions between doped hole carriers are very weak for 
two reasons: One is the low concentration of hole carriers in the underdoped regime 
and the other is the wave functions of hole carriers with up and down spins in a CuOe 
octahedron occupying aig and big orbitals, respectively, as seen in eq. (1) and Fig. 1. 
For these reasons, we have neglected the electron-electron interactions between doped 
holes in the effective Hamiltonian eq. (2). 

By replacing the localized spins SiS in iJcx with their average {S) in the mean-field 
sense, we can calculate the change in the total energy upon moving a hole from an 
aig orbital state in Hund's coupling spin triplet at Cu site i to an empty big orbital 
state in the Zhang- Rice spin singlet at the neighboring Cu site j. In the first step, the 
hole moves from Cu site i to infinity. The change in the total energy in the mean field 
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approximation is equal to ^aj^ + i-^'a*g- In the second step, the hole moves from infinity 

to an empty big orbital state at Cu site j to form the Zhang-Rice singlet. The change 

in the total energy in the second step is equal to ^big — f-^big- As a result, the change 

in the total energy by the transfer of the hole from the occupied a^g orbital state at Cu 

site i to the empty big orbital state at Cu site j is 

Here, ej*^ and e^_^ represent the effective one-electron energies of the aig and big 
orbital states including the exchange interaction term H^x, respectively. Thus, the en- 
ergy difference {e^J^ — £^big ) corresponds to the energy difference between the aig and 
big floors in the second story in Fig. 2. 



transfer interaction 




Fig. 3. (Color online) Simple explanation for coexistence of metallic state and local AF order in two- 
component theory (a*g and big orbital states), by taking neighboring Cu sites i to j as an example, 
where localized spins are in antibonding b^^g orbitals. Schematic pictures of Hund's coupling triplet 
^Big and Zhang-Rice singlet ^Aig are also shown at the bottom. The energy in this figure is taken as 
the hole energy. 

Now, let us estimate the energy difference ejf^ - e^_^ using the values of the 
parameters in the effective Hamiltonian eq. (2). The values of the parameters in 
Hamiltonian eq. (2) have been determined in the case of LSCO in ref. 6 (see also 
ref. 7). They are J = 0.1, K,*^ = -2.0, Ky,,^ = 4.0, ta-^a^^ = 0.2, tb.gbig = 0.4, 
tajgbig = A/^ajga*gibigbig ~ 0.28, ^a'^ = 0, and ^big = 2.6 in units of eV, where K^*^ 
and -ft'big are taken from first-principles cluster calculations for a CuOe octahedron in 
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LSCO/"'^"-' and the tmn are obtained by band structure calculation."'^-' The difference 

in one-electron energy between the a^g and big orbital states in a CuOe octahedron for 

a certain x has been determined so as to reproduce the difference in the lowest state 

energy between Hund's coupling spin-triplet state and the Zhang-Rice spin-singlet state 

for the same x in LSCO calculated by Multi- Configuration Self-Consistent Field (MC- 

SCF) cluster calculations which include the anti-JT effect.^) 

Thus, the calculated e^t^ - e-^^^ is 0.1 eV in the case of the optimum doping {x = 

0.15). Then, by introducing the transfer interaction of ta*gbig = 0.28 eV, a coherent 

metallic state in the normal phase is obtained in the presence of the local AF order for 

the underdoped regime. This situation is schematically shown in Fig. 3. 

2.3 Features of the many-body effect including energy bands and Fermi surfaces of 
underdoped LSCO coexisting with the AF order 

In the previous subsection, we have shown that the effective Hamiltonian eq. (2) 
for the K-S model can lead to a unique metallic state in the normal phase, which re- 
sults in the coexistence of a superconducting state and an AF order below Tc. In 1994, 
Kamimura and Ushio calculated the energy bands and Fermi surfaces of underdoped 
LSCO in the normal phase on the basis of the effective Hamiltonian eq. (2), by treating 
the fourth term H^x in the effective Hamiltonian eq. (2) by the mean-field approxima- 
tion, that is, by replacing the localized spins SiS with their average {S}.^'^'^ Thus, the 
effect of the localized spin system was dealt with as an effective magnetic field acting 
on hole carriers. As a result, Kamimura and Ushio separated the localized hole-spin 
system in the AF order and the hole carrier system from each other, and calculated 
the "one-electron type" energy band for a carrier system assuming a periodic AF or- 
der. Here, "one-electron type" means the inclusion of many-body effects in the energy 
bands. That is, the effect of the exchange interactions between carriers and localized 
spins is included in the sense of the mean field approximation. 

In Fig. 4, the calculated many-body effect including energy band structure for up- 
spin (or down-spin) doped holes in LSCO is shown for various values of the wave vector 

— * 

k and symmetry points in the AF Brillouin zone, where the AF Brillouin zone is adopted 

because of the coexistence of a metallic state and the AF order; it is shown on the left 
side of the figure. Here, note that the energy in this figure is taken as the electron 
energy but not as the hole energy. Furthermore, the Hubbard bands for localized hi 
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Fig. 4. Many-body effect including band structure*' for up-spin (or down-spin) dopant holes in 
undcrdoped LSCO above T^. The highest occupied band is marked by the #1 band (right) and the 
AF Brillouin zone (left). The A point corresponds to (7r/2a, 7r/2a, 0), while the Gi point corresponds 
to (7r/a,0, 0), at which a saddle-point singularity appears. 

holes, which contribute to the local AF order, are separated from this figure and do not 
appear in this figure. 

In undoped La2Cu04, all the energy bands in Fig. 4 are fully occupied by electrons 
so that La2Cu04 is an AF Mott insulator, consistent with experimental results. In 
this respect, the present effective energy band structure is completely different from 
the structure of ordinary LDA energy bands. When Sr is doped, holes begin to 
occupy the top of the highest band in Fig. 4 marked by #1 at the A point, which 
corresponds to (7r/2a, 7r/2a, 0) in the AF Brillouin zone. At the onset concentration of 
superconductivity, the Fermi level is located slightly below the top of the #1 band at 
A, which is slightly higher than the Gi point. Here, the Gi point in the AF Brillouin 
zone lies at (tt/c, 0, 0) and corresponds to a saddle point of the van Hove singularity. 

On the basis of the calculated band structure shown in Fig. 4, Kamimura and 
Ushio^'^^ calculated the Fermi surface (FS) for the underdoped regime of LSCO. The 
calculated FS in the underdoped regime is composed of four Fermi pockets of extremely 
fiat tubes. The projected two-dimensional (2D) picture of the four Fermi pockets around 
the A point, (7r/2a, 7r/2a) and the other three equivalent points in the momentum space 
is shown in the antiferromagnetic Brillouin zone in Fig. 5. The total volume of the four 
Fermi pockets is proportional to the concentration of the doped hole carriers. Thus, the 
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Fig. 5. (Color online) Fermi surface (FS) for up-spin (or down-spin) dopant holes in underdoped 
LSCO above Tc. The FS consists of four Fermi pockets around the A point, (7r/2a, 7r/2a, 0), and 
the other three equivalent points (the nodal region) in momentum space. The figure shows the two- 
dimensionally projected Fermi pockets. The point A represents one of the electronic states. 



feature of Fermi pockets constructed from the doped holes shown in Fig. 5 is consistent 
with Luttinger's theorem in the presence of AF order. 

In 1996 and 1997, respectively. Mason et al?^^ and Yamada et al?^^ independently 
reported the magnetic coherence effects on the metallic and superconducting states 
in underdoped LSCO, determined by neutron inelastic scattering measurements. Since 
then, a number of papers suggesting the coexistence of local AF order and superconduc- 
tivity in cuprates as a result of neutron and NMR experiments have been published. ^''"^^^ 

The Fermi surface structure in Fig. 5 is completely different from that of the single- 
component theory, in which the FS is large. Recently, Meng et al. have reported the 
existence of the Fermi pocket structure in the ARPES measurements of underdoped 
Bi2Sr2_^La^Cu06+5 (La-Bi2201).^*^'' Their results are clear experimental evidence of 
our Fermi pocket structure for underdoped LSCO predicted in 1994.*^ 

In 1997, Anisimov et al. calculated the energy band structure of the ordered alloy 
La2Lio.5Cuo.5O4 by the LDA+U method, '^^^ and they showed that a fairly modest reduc- 
tion in the apical Cu-0 bond length is sufficient to stabilize Hund's coupling spin triplet 
state with dopant holes in both big and a*g orbitals. Their calculated result supports 
the K-S model. 
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3. Cakulaliuii uf ARPES SptjcLia Based on llie K-S Model and llie Cun- 

clusion of the Absence of a Pseudogap 

Recently, considerable attention has been paid to the phenomenological idea of the 
pseudogap. When a portion of the Fermi surface in cuprates was not observed in the 
ARPES experiments, the idea of the pseudogap was proposed as a type of gap for 
truncating the FS in a single-particle spectrum. ^^'^^^ The disconnected segments of the 
FS are called the "Fermi arc".^^'^^'^^^ Further ARPES experiments showed that such 
a pseudogap develops below a temperature denoted T*, which depends on the hole 
concentration x in the underdoped regime of cuprates; thus, we write T*{x) hereafter. 
T*[x) decreases with increasing hole concentration x and disappears at a certain con- 
centration Xo in the overdoped region. In this section, on the basis of the K-S model, 
we clarify the origins of the pseudogap and T*{x). 

3.1 Calculation of the photoemission intensity and clarification of the origin of the 
observed two distinct gaps 

Below Tc, the hole carriers in the Fermi pockets shown in Fig.5 form Cooper pairs, 
contributing to the formation of a superconducting state, and a superconducting gap 
appears across the Fermi level. This feature is consistent with Uemura et a/.'s plot.^^^ 
In Fig. 6(a), the d-wave node below predicted by the K-S modeP°'^^^ is schematically 
shown as dots, and the d-wave superconducting density of states is schematically shown 
in Fig. 6(b). Here, note that the AF order still coexists with a superconducting state 
below Tc so that we can use the same AF Brillouin zone, as shown in Fig. 6. 

On the other hand, in the antinodal region, the states occupied by electrons that do 
not participate in the formation of superconductivity still exist below T^.. As an example 
of such states, the state A is shown in Fig. 6(a), and the state corresponding to A above 
Tc is also shown in Fig. 5. Then, real transitions of electrons from the occupied states, 
say, the state A, below the Fermi level in the #1 energy band in Fig. 4 to a free- 
electron state above the vacuum level occur by photoexcitation both above and below 
Tc around the Gi point (tt/c, 0,0) and other equivalent points in momentum space. 
These transitions appear in the antinodal region in momentum space. 

Such a transition is shown in Fig. 7.^^^ Let us consider the case in which an electron 
in the occupied state A with energy Si and momentum ki in the ^1 energy band below 
e-p is excited to a frcc-clcctron state with the bottom of the energy dispersion at {ep—eo) 
in a crystal by a photon with energy hu, where we use the suffix i to emphasize the 
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Fig. 6. (Color online) Sketch of d-wave superconductivity in K-S model below Tc- (a) Change of 
Fermi pockets in nodal region to d-wave nodes below T^. (b) d-Wave superconducting density of 
states. 



initial state of the transition in the crystal. Thus, the final state of the transition with 
energy in the crystal is expressed as 

Si = —{kl + k]) + sy - So, (4) 

where k\\ and k_[_ are the momenta of the photoexcited electron parallel and perpendic- 
ular to the crystal surface, respectively. As shown in Fig. 7, the energy conservation for 
this excitation process from the initial state \i) to the final state |/) in the crystal is 
expressed as hv = Sf — £\ for a photon of energy hv. 

When an electron is ejected into the vacuum level of the crystal by a photon with 
energy hv, it acquires kinetic energy. Through ARPES experiments, we measure the 
kinetic energy of photoelectrons emitted in vacuum. We define the kinetic energy of 
such photoelectrons emitted in vacuum as ^kin, where £kin = (/?.^/2m)(i^'| + K^). 

On the other hand, the binding energy of the electron in the initial state, E-q, is 
introduced as a new variable instead of £{ki). E-q is defined as 

Eb = -e{h) + ep. (5) 

The following equation also holds for ^f: 

£i = Skin + W + eF, (6) 
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Fig. 7. (Color online) Energetics of photoemission process. 



where W is the work function of the crystal (see Fig. 7). 

In the ARPES experiment, when a photoelectron is emitted from a crystal in vac- 
uum through a surface, it is assumed that the momentum parallel to the surface is 
conserved: fcy = K\\. Now, note that the #1 band in Fig. 4 has been calculated by the 
mean field approximation for the fourth term in the Hamiltonian (2). The important 
consequence of this approximation is that, having taken into account the strong spin 
exchange interaction in the mean field approximation, the probability of removing an 
electron in the state with momentum k\ and energy Si in the t^I energy band to the free- 
electron state in vacuum can be treated in a framework similar to that for single-particle 
photoexcitation. 

As a result, the EDCs in the ARPES experiments corresponding to the transition 
from the occupied states in the many-body effect including energy band in Fig. 4 to the 
free-electron band can be calculated using the following formula for the photoemission 
intensity I{k,u): 



Furthermore, pi{u) is the density of final states and |Mg| is the squared one-electron 




Here, A{k, u) is the spectral function that gives the probability of removing or adding an 
electron at {k, u), where u is the electron energy relative to the Fermi level. It is related 
to the imaginary part of the one-electron Green's function; A{k, u) = —{l/iT)lmG{k.u). 



2 
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transition matrix element. It is clear from b'ig. 7 that A[k,uj) gives the highest prob- 

ability when huj is equal to {hu — £kin — W), where hu — e{ki) - e{ki). 

By taking account of the hfetime effects due to the finite size of a metallic state, the 

deviation from the mean field approximation, and other factors, the spectral function 

A{k,uj) is given by 

A{k,u;) = (I/tt)- \ (8) 

where 5 denotes the lifetime effects, and the momentum dependence in Si is expressed 
as e{ki) explicitly. 

The density of final states Pf(£:km) in the EDCs is defined from the dispersion of 
the momentum of a photoexcited electron perpendicular to the crystal surface in the 
crystal, k±, as, 

dfcx 



Pf(£kin) 



(9) 



d^kin 

Using eqs. (4) and (6) with the conservation of momentum of a photoexcited electron 
parallel to the crystal surface, k\\ — Kh , the density of final states is obtained as. 



Pf(£kin) = , (10) 

2 J{hy2m) (skin + V- {hy2m)klj 



where V = W + So is the inner potential. This result agrees with the result derived by 
Mizokawa.^^-' 

Since much of the ARPES EDC data is expressed as a function of the binding energy 
^^B, we express eqs. (8) and (10) in terms of Eb- For this purpose, we first insert eq. (5) 
into eq. (8), and simultaneously replace hut in eq. (8) by {hv — £kin — W). As a result, 
eq. (8) can be written as. 

Furthermore, using the expression for the inner potential, V — W + Eq, and the 
energy conservation relation in Fig. 7 given as 

hiy = £kin + W + Es, (12) 
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eq. (10) can be expressed as 

Pf(£kin) = — , ^ =, (13) 

where kq is the component of ki paraUel to the crystal surface. 

Using eqs. (7), (12), and (13), we have calculated the photoemission intensity I{k, cu) 

— * 

as a function of E-Qih). In performing the numerical calculations, we have considered 
that the photon energy (hi/) range in synchrotron radiation experiments is 10 to 100 
eV and the kinetic energy range of the photoelectron is also 10 to 100 eV.^^^ Since the 
width of the energy dispersion of the #1 energy band in Fig. 4(a) is about 1 eV, we 
notice that the E-q range is up to 1 eV. For S, whose inverse gives a measure of the 
lifetime broadening in the #1 band, we assume 100 meV on the basis of the discussion 
in the subsequent section. 

In this context, we choose 15 eV for hu, 10 eV for £kin, 3.5 eV for W,"^^^ and 4.5 eV 
for the inner potential V ~ So + W in the present numerical calculations. As regards 
h^k1^^/2m, we choose the center of the antinodal region, i.e., the Gi point or the edge 
of the AF Brillouin zone in Fig. 4, for the ith point, because the antinodal region is 
narrow around the Gi point, so that h^k^^/2m does not change much upon varying the 
ith point. By adopting the empty lattice test for the free-electron energy bands, we 
estimate h^k1^^/2m to be 3 eV for i — Gi. 

— * 

The calculated I{k, cu) with the values of the above parameters is shown as a function 
of Eb in Fig. 8(a). Since etm is very large, a divergent point in the density of final states 
Pf{uj) appears at a large E^. Thus, the photoemission intensity I{k,Lj) shows a feature 
of a broad hump, reflecting a peak in the spectral function A{k,u) given by eq. (11), as 
seen in Fig. 8(a). This trend is consistent with the experimental results of the ARPES 
spectra of underdoped Bi2212 samples below Tc in the antinodal region by Tanaka et 
al.,^^^ although the shape of the broad hump is shghtly different. 

From the ARPES spectra in the nodal region shown in Fig. 6(b), which was predicted 
from the d-wave superconductivity due to the K-S model, and those in the antinodal 
region shown in Fig. 8(a), we can conclude that the features of ARPES spectra below 
Tc are theoretically as follows: ARPES spectra consist of a coherent peak due to the 
superconducting density of states that appears in the nodal region around the A point 
and a broad hump that appears in the antinodal region, which corresponds to the 
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Fig. 8. (Color online) Calculated ARPES spectra of LSCO and their comparison with experimental 
results of Bi2212. (a) Calculated photoemission intensity as a function of binding energy Eb{— E—Ep). 
(b) Calculated energy difference |e(Gi) — eF(2^)| (antinodal transition energy) as a function of hole 
concentration a; at T = OK. Experimental results of Tanaka et al}^^ are shown by dots. 

real transition of electrons from the occupied states below the Fermi level ep to a 
free-electron state above the vacuum level. These theoretical results of ARPES EDCs 
are similar to the experimental ones reported by Tanaka et al. for Bi2212/^) where 
the experimental results revealed two distinct energy gaps in the nodal and antinodal 
regions exhibiting different doping dependences. Thus, we designate the broad hump 
in the antinodal region as an "antinodal transition", where £f varies with the hole 
concentration a;, so that we write e-p{x) hereafter. 

Furthermore, Tanaka et al. reported the doping dependence of the position of the 
hump, which is determined from the second derivative of the spectra in the antinodal 
region around the Gi point for the three underdoped samples in Bi2212. We compare 
this experimental result of Bi2212 with the calculated doping dependence of the bind- 
ing energies Eb at the Gi point for LSCO, which correspond to the energy difference 
|£:(Gi) — e-p{x)\. We call this energy difference the "antinodal transition energy". 

Since the shape of the density of states (DOS) for the highest conduction band does 
not depend on the type of cuprate material, we can compare in Fig. 8(b) the calculated 
doping dependence of the antinodal transition energy for LSCO (solid lines) with the 
experimental results of the antinodal gap of Bi2212 in ref.l6, which are shown as dots in 
the figure. As seen in Fig. 8(b), the agreement between the theory and the experiment 
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IS remarkably guuu. Fruiii Ihis quunLiLaLivt! Hgr(j(jiiieiiL, we can cuiiclude LliaL amuug Lne 

observed two gaps below Tc, the gap associated with the antinodal regime corresponds 
to the real transitions of electrons from the occupied states below the Fermi level to a 
free-electron state above the vacuum level, while the other gap associated with the near- 
nodal regime corresponds to the superconducting gap created on Fermi pockets. From 
the excellent agreement between the present theoretical results and the experimental 
results of Tanaka et al, we can conclude that the real transitions of electrons from 
occupied states below the Fermi level to a free-electron state above the vacuum level 
by photoexcitation appear in the antinodal region in underdoped cuprates so that the 
introduction of a pseudogap is not necessary. 

Recently, Yang et al.^^^ have suggested from their ARPES experiments on Bi2212 
that the opening of a symmetric gap related to superconductivity occurs only in the 
antinodal region and that the pseudogap reflects the formation of preformed pairs, 
in contrast to the ARPES experimental results reported by Tanaka et al}^^ In the 
present theory, we have clearly shown in Figs. 6 and 8 that, in the ARPES experiments, 
a peak related to superconductivity appears only in the nodal region and that the 
spectra in the antinodal region correspond to photoexcitations from occupied states 
below 5f to a frcc-electron state above the vacuum level. If the antinodal region in 
ref. 44 is the region around the Gi point in the present paper, their suggestion is in 
disagreement with our theoretical results. Finally, we should remark that any proposed 
theory must explain both the doping and temperature dependences of ARPES spectra 
in the underdoped regime consistently. From this standpoint, we will investigate the 
temperature dependence of the calculated broad hump in ARPES EDCs theoretically 
in the next subsection. 

3.2 Physical meaning ofT*{x) and the temperature dependence of ARPES spectra 

To calculate the temperature dependence of the antinodal transition energy, first we 
would hke to clarify the physical meaning of T* (x) . When a hole concentration x is fixed 
at a certain value in the underdoped region and the temperature increases beyond Tf., 
in the normal phase, the local AF order constructed by superexchange interaction in a 
Cu02 plane is destroyed by thermal agitation, and thus a phase showing the coexistence 
of a metallic state with the Fermi pockets and the local AF order diminishes gradually. 
As a result, an electronic phase consisting of a large FS without the AF order is mixed 
with a phase of the K-S model. Finally, at a certain temperature, a uniform phase 
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appear in the underdoped regime. This temperature is defined as T*{x). Thus, the 

phase of the Fermi pockets coexisting with the local AF order in the K-S model holds 

only below T*{x). We designate the phase of the Fermi pockets in the K-S model as 

the "small FS" phase and the electronic phase consisting of a large FS without the AF 

order as the "large FS" phase. Hereafter, the former and latter are abbreviated as the 

SF and LF phases, respectively. In this context, one may consider that a phase below 

T*{x) is a mixed phase of the SF and LF phases in the underdoped regime; thus, T*{x) 

represents a crossover from the mixed phase to the LF phase. 

To calculate T*{x) on the basis of the K-S model, one must take account of the effect 
of thermal agitation in the system of Cu localized spins in the AF order (the first story 
in Fig. 2). However, such calculation is possible only for a finite system, as Hamada and 
coworkers have shown.^^'^^^ In this context, we calculate T*{x) approximately, neglecting 
the effect of thermal agitation in the system of Cu locahzed spins. 

For this purpose, let us introduce a quantity that defines the difference between the 
free energies of pure LF and SF phases: 

AF{T,x) = Fi^f{T,x)-Fsf{T,x), (14) 

where Fi^f{T,x) and Fsf{T,x) are the free energies of the LF and SF phases, respec- 
tively. Here the free energy F(T, x) is defined as 

F{T,x) = E{T,x) -TS{T,x), (15) 

where E{T,x) and S{T,x) are the internal energy and entropy of each phase, respec- 
tively. These quantities are calculated from 

/oo 
£p(£)/(£,/x(x))d£, (16) 
-oo 

and 

/oo 
/(£,Ma:))ln/(£,Mx)) 
-oo 

+ {l-/(£,Ma;))}ln{l-/(£,M^))} 

p{e)de, (17) 

where jji{x) is the chemical potential of each phase, p{e) is the DOS for each phase, and 
f{e,p{x)) is the Fermi distribution function at energy s and chemical potential p{x). 
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AF{f*{x),x) = 0. (18) 
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Fig. 9. (Color online) (a) Calculated temperature dependence of antinodal transition energy for x 
= 0.05, 0.1, and 0.15. Experimental data obtained by Norman et al.^'^^ (triangles) and Lee et al.^^^ 
(squares) are shown, (b) Calculated result of T*{x) as a function of the hole concentration x 



Kamimura et al. calculated the electronic entropies for the SF and LF phases of 
LSCO.^^ According to their results, the difference in electronic entropy between the 
SF and LF phases increases with increasing hole concentration x in the underdoped 
regime. Using this result, we have calculated T*{x) from eq. (18) as a function of x, 
instead of T*(x). In doing so, we have introduced two parameters, T*(x = 0.05) and 
Xo, where Xo is the critical concentration that satisfies T*(xo)=0. Here, T*{x = 0.05) 
represents a quantity related to the energy difference between the phase of the doped AF 
insulator and the LF phase at the onset concentration of the metal- insulator transition 
{x = 0.05). T*{x = 0.05) is chosen to be 300K. On the other hand, the physical meaning 
of T*{xo) = can be explained as follows: When the hole concentration exceeds the 
optimum doping level {x = 0.15) for LSCO and enters a slightly overdoped region, the 
local AF order via superexchange interaction in a Cu02 plane is destroyed by an excess 
of hole carriers. Thus, the K-S model does not hold at a certain concentration Xq in the 
overdoped region, and hence the small FS in the K-S model changes to a large FS. Thus, 
T*{x) vanishes at Xo- From the analysis of various experimental results, we choose Xo 
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= U.iV for LSCO>^'^°) 'i'he calculated result of T*{x) with T*{x = U.Ub) = 3UUK and 

Xo — 0.17 is shown as a function of x in Fig. 9b. 

Prom the present result, we can say that the area below T*{x) in the underdoped 
regime represents the region in which the normal (metallic) phase above and the su- 
perconducting phase below Tc coexist with the local AF order. In a real system, a region 
of a mixed phase consisting of the SF and LF phases appears between T*{x) and T*{x) 
owing to the dynamical interaction of the fourth term in the effective Hamiltonian (2). 
Thus, T*(x) always appears above T*(a;).^^-' 

Under this circumstance, it is clear that the antinodal transition energy defined 
by |£(Gi) — eF{x)\ appears at temperatures below T*{x) and vanishes at T*{x). By 
using T* {x) instead of T* {x) , we calculate the temperature dependence of the antinodal 
transition energy using eqs. (14)-(17). The calculated results for three concentrations, 
i.e., X — 0.05, 0.10, and 0.15, of LSCO in the underdoped-to-optimallydoped region are 
shown in Fig. 9(a) as functions of temperature, where T*(0.05) = 300K is used. As seen 
in the figure, the antinodal transition energy increases slightly with temperature up 
to T*{x) and vanishes suddenly at T*(x). These calculated results are compared with 
the experimental results of the underdoped sample of Bi2212 in refs. 47 and 48, which 
are indicated by triangles and squares, respectively, in Fig. 9(a). As seen in the figure, 
the agreement between the theory and the experiment is remarkably good, although 
the present theoretical treatment is not rigorous in the sense that, in the calculation of 
Fsf{T,x), the dynamical interplay of a metallic state (the second story in Fig. 2) and 
Cu locahzed spins in the AF order (the first story) is not taken into account. 

4. Spatially Inhomogeneous Distribution of Fermi Pocket States and Large 
Fermi Surface States due to the Finite Size Effects 

4-1 Finite size effects of metallic state on Fermi surface 

According to the results of neutron inelastic scattering experiments by Mason et 
al.^^^ and Yamada et al.,'^^^ the AF spin-correlation length in the underdoped region 
of LSCO is finite. In the underdoped regime of LSCO, it increases as the Sr concentration 
increases from x = 0.05 in LSCO, the onset of superconductivity, and reaches about 
50A or more at the optimum doping level {x = 0.15). In this subsection, we discuss the 
effects of the finite size of the AF spin-correlation length on the structure of the Fermi 
pockets shown in Fig. 5. According to the K-S model in Fig. 2, in the spin-correlated 
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regiuli a aup(jd liul(j iii Llie uiiderduped regime uf LSCO can lUiieruLt; cuhereiiLly by 
taking the a^g and big orbitals alternately in the presence of the local AF order without 
destroying the AF order. 

In the case of a finite spin-correlated region, one may think that there are frustrated 
spins at the boundary between the spin-correlated region of the AF order and the 
region of the "resonating valence bond" (RVB) state proposed by Anderson^) without 
hole carriers. Here, the frustrated spins mean that the localized spins at the boundary 
are not in the AF order, but directed parallel to each other. Suppose that one of the 
frustrated spins in a parallel direction at the boundary changes its direction from parallel 
to antiparallel by the fluctuation effect in the 2D Heisenberg AF spin system during 
the time of Tg defined by = h/J, where J is the superexchange interaction (~0.1 
eV). At the time of r^, on the other hand, hole carriers at the Fermi level can move 
with the Fermi velocity inside the spin-correlated region of the AF order. The traveling 
time of a doped hole at the Fermi level over an area of the spin-correlation length is 
given by rp = Xs/vf, where vp is the Fermi velocity of a doped hole at the Fermi level. 
In the case of underdoped LSCO, is 6 x 10~^^ s. Since vp is estimated to be 2.4 
X 10^ m/s from the dispersion of the #1 band in Fig. 4, rp is 2 x 10~^^ s for the 
underdoped region of x = 0.10 to x = 0.15 in LSCO, where for the spin-correlation 
length Xg Sit X = 0.15, wc have chosen 50 A. Thus, rp becomes much longer than r^. As 
a result, the frustrated spins on the boundary change their directions from parallel to 
antiparallel before a hole carrier in the spin-correlated region of the AF order reaches 
the boundary. Thus, a metallic state for a doped hole becomes much wider than the 
observed spin-correlation length by the passing of a doped hole through the boundary 
without spin scattering. In this way, a metallic state is surrounded by RVB states and 
the spatial distribution of the metallic states is inhomogeneous. 

Finally, we explain why we have chosen 100 meV for S in calculating the photoe- 
mission intensity shown in Fig. 8(a). For example, the initial state of photoexcitation 
in ARPES near the Gi point is either a component of Hund's coupling triplet ^Big or 
Zhang-Rice singlet ^Aig shown in Fig. 3 in the ^1 band.^-* Thus, if the local AF order 
between neighboring Cu sites in Fig. 3 is destroyed, the calculated result in Fig. 8(a) 
may not be valid. This is the reason why S is the same order of magnitude as the inverse 
of Tg, that is, the superexchange interaction J (~0.1 eV). 
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'^rr2 — Onyvii oj iht coextstence uj a local AF order and a iiietalhc state: Tlit kmettc- 

energy-driven mechanism 

Concerning the finite system of cuprates, Hamada et al.'^^^ and Kamimura and 
Hamada^^^ attempted to determine the ground state of the effective Hamiltonian (2) 
for the K-S model by carrying out the exact diagonahzation of the Hamiltonian (2) 
using the Lanczos method for a 2D square lattice system with 16 (4 x 4) locahzed 
spins with one and two doped holes, respectively. As a result they clarified that, in the 
presence of hole carriers, the localized spins in a spin-correlated region tend to form an 
AF order rather than a random spin-singlet state, and thus hole carriers can lower the 
kinetic energy by itinerating in the lattice of the AF order (the first story in Fig. 2). 
This is the mechanism leading to the coexistence of a metallic state and a local AF 
order in the K-S model. 

Generally, a hole-carrier in the spin-correlated region of the AF order can propagate 
through the boundary of the spin-correlated region with the above-mentioned mecha- 
nism of the K-S model; hence, the region of a metaUic state coexisting with the AF 
order becomes much wider than the observed spin-correlated region. In fact, Kamimura 
et al. estimated the length of the metaUic region at the optimum doping level of LSCO 
to be about 300 A from the at the optimum doping level. Recently, an idea similar 
to ours with regard to the decrease in the kinetic energy has been proposed by Wrobel 
and coworkers, who have shown that the decrease in the kinetic energy is the driving 
mechanism that induces superconductivity.^^' 

5. Remark on a Phase Diagram for Underdoped Cuprates 

From the calculated results shown in Fig. 9(b), we would like to comment on the 
T \s X phase diagram for cuprates shown in Fig. 10, for which it has been said that a 
pseudogap state exists below the temperature T*{x). According to our calculations in 
previous sections, the SF phase constructed from Fermi pockets appears in the pres- 
ence of the local AF order below T*{x) in the underdoped region. However, when the 
temperature increases at a fixed concentration in the underdoped regime, the AF order 
is destroyed gradually with increasing temperature, and thus the K-S model does not 
hold slightly below T*{x). On the other hand, when the hole concentration increases 
at a fixed temperature, the AF order is destroyed by overdoped holes. Thus, the K-S 
model does not hold at a certain hole concentration. The thermal effect and excess hole 
density effect cause a mixing of the SF and LF phases, as explained in §4. 
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In Lliis conLexL, we would like Lo poinL OLiL LliaL the area below T*{x) and above 

in Fig. 10 is not a pure phase but an inhomogeneous distribution of Fermi pocket states 
and large FS states. Thus, we can say that, when the temperature increases from T^. 
at a fixed hole concentration in the underdoped region, the population showing Fermi 
pockets decreases while that showing large FS increases with increasing temperature. 
Finally, when the temperature exceeds T*{x), a uniform LF phase appears as a metallic 
state. Thus, we may say that T*{x) is a crossover line from the inhomogeneous mixed 
phase of Fermi pockets and large FS to the LF phase. 

This result can explain the strange temperature evolution of a Fermi arc observed 
by Norman et al.^^^ and Kanigel et al.^"^^ Furthermore, in the superconducting phase 
below Tc and T*{x) in Fig. 10, an s-wave component of superconductivity originating 
from the LF phase may be mixed with the d-wave superconductivity. Such a mixing 
effect was experimentally reported by Miiller.^^^ In this context, we conclude that T*(x) 
represents a crossover from the SF phase to the LF phase rather than a phase boundary 
between the pseudogap phase and a metal. 



UK) 




Mainly d-wave superconducting phase 



Fig. 10. (Color online) New explanation for phase diagram of LSCO. Since T*{x), which is higher 
than T*(a;), has not been calculated, its numerical values are not shown on the vertical axis. 

Furthermore we can predict that the spin susceptibility will show 2D-like AF features 
mainly below T*{x) and Pauli-like temperature-dependent behavior above T*(x). We 
find that this prediction is also consistent with the experimental results for LSCO.^^'^*^^ 
In this context, it should be emphasized that the K-S model is shown to explain suc- 
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cessfuUy not only the ARPES experimental results^^'^'"'^''''^"-* but also a number of other 

experimental results such as NMR results showing the coexistence of a superconduct- 
ing state and AF order, polarized X-ray absorption spectra,^^'^^^ site-specific X-ray 
absorption spectroscopy,^^^ anomalous electronic entropy/' and d-wave superconduc- 
tivity. Theoretically, the K-S model is also supported by LDA -|- U band calcula- 
tions,^^) as already mentioned in §2.5. 

6. Conclusions and Concluding Remarks 

In this study we have shown on the basis of the K-S model how the interplay of Mott 
physics and JT physics plays an important role in determining the superconducting state 
as well as the metallic state of underdoped cuprates. In connection with the interplay 
of JT physics and Mott physics, the following important results have been obtained 
in this study: It has been clarified on the basis of the K-S model that the concept 
of the pseudogap discussed theoretically^^' and reported by ARPES, STM, and 
tunneling experiments below T*{x) in underdoped cuprates^^'^^'^^^ can be explained 
by the occurrence of Fermi pockets in the underdoped region without the pseudogap 
hypothesis. We have shown that the appearance of the broad hump observed in the 
antinodal region in ARPES can be explained by the real transitions of photoexcited 
electrons from the occupied states in the highest conduction band in the antinodal 
region to a free-electron state above the vacuum level. Furthermore, we have shown that 
the physical meaning of T* represents a crossover line from an inhomogeneous phase 
consisting of Fermi pockets and large FS in the normal state to a phase consisting of a 
large FS. 

Finally, several remarks are made on the small FS and shadow bands in the un- 
derdoped regime of cuprates. In 1996, Wen and Lee developed a slave-boson theory 
for the t-J model at finite doping level, and showed that Fermi pockets at low doping 
continuously evolve into a large FS at high doping concentrations.^^) Although their 
theoretical model is different from the K-S model, it is interesting to find that they 
obtained a result similar to the result predicted using the K-S model in 1994 with re- 
gard to the change from a small FS to a large FS with increasing hole concentration. 
Recently, a proposal was made to reconcile the experimental result of the coexistence 
of antiferromagnetism and super conduct ivity.^^-* Furthermore, in relation to the small 
FS, the idea of a shadow FS was proposed as a replica of the main FS transferred using 
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Q = [71 /a, 71/ a) by Kampf and Schrieffer theoretically^'''^ and by Aebi et al. experi- 

mentally.^^^ Investigating the vahdity of the idea of the shadow FS experimentally, the 

observation of shadow bands in ARPES spectra has been reported. ^^^^^ Responding 

to the problems of the shadow FS and shadow bands from the standpoint of the K-S 

model, it should be emphasized that Fermi pockets in the metallic state calculated from 

the K-S model have been derived as a result of the interplay of JT physics and Mott 

physics; thus, the origin of Fermi pockets is different from that in a single-component 

theory. Therefore, the Fermi pockets shown in Fig. 5 are neither the shadow FS nor 

related to the shadow bands. 
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